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Pursuing the original idea proposed in our previous paper (Paper I), we improve the 
method to determine the shape of the initial curvature perturbation spectrum P(k) from 
the CMB data. The thickness of the last scattering surface (LSS) and the integrated 
Sachs- Wolfe (ISW) effect, which we neglect in Paper I, are taken into account and an it- 
erative method is newly developed. The new method can reproduce the primordial power 
spectra with a high accuracy, given the correct values of the cosmological parameters. 
Conversely, there appear spurious peaks and dips in the reconstructed power spectrum 
if we use the cosmological parameters slightly different from the true values, while there 
appear regions of negative P(k) in some cases if we use substantially different values. In 
other words, the tacit assumption that the cosmological parameters can be determined 
for an assumed initial spectrum is verified by our reconstruction method. In addition, it 
could be a new tool to constrain the cosmological parameters without recourse to models 
of the primordial power spectra. 

PACS: 95.30.-k; 98.80.Es 



I. INTRODUCTION 

Cosmic microwave background (CMB) is now playing a central role in the precision cosmology and 
enables us to extract wealth of information on the cosmological parameters and the primordial universe. 
The past CMB experiments suggest that our universe is consistent with spatially flat ACDM universe 
with scale- invariant initial spectrum [1], which is predicted by conventional slow-roll inflation models [2]. 
Nevertheless these results are not too restrictive and other cosmological models are still possible under 
the present accuracy of the observational data. In the near future, we will attain unprecedentedly precise 
data by MAP [3] and Planck [4] experiments. Then, we will be able to discuss the results much more 
quantitatively than we can now. 

In particular, on the analysis of the CMB anisotropy, it is more desirable to deal with the initial 
spectrum of perturbation as an arbitrary function, although most of the past analyses of observational 
data assumed a power-law spectrum. On the theoretical side, a variety of generation mechanisms of 
broken scale- invariant spectra have been proposed, even in the context of inflationary cosmology [5]. 
For this reason or others, there have been some attempts to reconstruct the initial spectrum from the 
data [6-8]. But some of them assumed that the initial spectrum is a piecewise power-law function or 
parametrised as a broken scale-invariant function [6], while others took only the Sachs- Wolfe (SW) effect 
[9] into account [7]. 

In our previous paper (hereafter Paper I) [10], we built up a basic framework to reconstruct the pri- 
mordial spectrum of curvature perturbations, P{k), out of CMB anisotropy data by solving a differential 
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equation for P(k), in which not only the SW effect but also the Doppler effect are taken into account. In 
Paper I, however, we considered only a simplified and rather unrealistic cosmological model for illustra- 
tion without any degrees of freedom of cosmological parameters. Furthermore, we neglected the thickness 
of the last scattering surface (LSS), which caused considerable errors in the angular power spectrum C^. 

In this paper, we improve the method to a more accurate one which can be applied to realistic models. 
The width of the LSS and Integrated Sachs- Wolfe term are newly included. We also investigate the 
dependence on the cosmological parameters. Our inversion method assumes that the true cosmological 
parameters are known. In reality, this is not the case. In fact, the standard argument is to use the CMB 
data to determine the cosmological parameters, with some assumptions on the initial power spectrum. 
However, we cannot apply this argument to our case, since we treat the initial spectrum as an arbitrary 
function. Thus, it is important to study the effect on the reconstructed spectrum as we vary the cosmo- 
logical parameters. Since our main purpose here is to complete solving the inverse problem for realistic 
models, we assume an ideal situation where the observational error does not exist. 

This paper is organized as follows. In Sec. II, we describe our inversion method in some detail, which 
is an improved version of the one presented in Paper I. In Sec. Ill, we apply our method to some model 
spectra, assuming the true cosmological parameters are known. We find that the new method can 
reproduce the primordial spectrum to quite accurately. In Sec. IV, we investigate the dependence of the 
reconstructed spectrum on the cosmological parameters. We find the cosmological parameters can be 
constrained severely without any assumption on the primordial spectrum. We discuss the implications 
of our method in Sec. V. 



II. INVERSION FORMULA 



First we write down basic equations. We use the same notation as Paper I. Working in the Fourier 
space and denoting the direction cosine between the wavenumber vector k and the momentum vector of 
the photon by /i, an integral form of the Boltzmann equation for the temperature anisotropy, 6(77, n, k), 
is given in the Newton gauge by [11] 

6(770, + *(%,*0 = J o {(eo + *- i /iV 6 )V(r ? ) + (*-$) e - T ( I ')} e 4fc ^-" )dr?. (1) 

Here, an overdot denotes a derivative with respect to the conformal time 77, and 770 is its present value, \& 
and $ are the Newtonian potential and spatial curvature perturbation on the Newton slices, respectively 
[12], and 60 is the monopole component of the multipole expansion, 

e( v ,k,v) = J2H) e Qi(v,k)p e (ri. (2) 
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The function V(r?) is called the visibility function and r(rj) is the opacity, which are given by 

Vfa) = T(rj)e- T ^ , rfa) = / ttf)drf, f = ax e n e a T , (3) 

where a is the cosmic scale factor, x e n e is the free electron density, and ot is the Thomson cross section. 
Note that we have neglected a term due to anisotropic stress in the integrand of Eq. (1). As we will see 
below, however, this causes no problem in our inversion method. Also, we note that W = — $ apart from 
the effect of small anisotropic stress due to photons and neutrinos. 

Under the thin LSS approximation adopted in Paper I, these functions are approximated by 

V(v) = 6(r,-*), e- T(r,) =%-^), (4) 

respectively [11], with 77* being the decoupling epoch when the visibility function is maximum. Hence, if 
we neglect the ISW term, 0£ (770 , k) is given by 

ef in (77o,fc) = {Q + ^)(r 1 ,,k){2£+l)j e {kd) + Q 1 (r],,k)(2£+l)j' e (kd), (5) 

where d = 770 — 77* is the conformal distance to the LSS. Assuming primordial fluctuation is adiabatic, 
both (60 + 'J) (t/*, k) and 61(77*, k) can be expressed as 

(60 + 10(77*, A;) =/(77*,fc)$(O,A0, 

6i(77*,A;)= 5 (77*,A ; )$(0,fc), (6) 
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where /(f?*, k) and g(r]*, k) are the transfer functions for the respective quantities. 

In Paper I, using these transfer functions, we derived a formula that relates the initial curvature 
perturbation spectrum and the CMB angular correlation function in a flat CDM universe under the thin 
LSS approximation: 

C(r)=3rC(r)+r 2 C'(r) = — / dkP(k) {f 2 (k)k 2 r cos kr + [2f 2 {k) + g 2 {k)] ksmkr} , (7) 

2^ Jo 

where P(k) = (|$(0, fc)| 2 ) is the initial spectrum, C(r) is the CMB angular correlation function, and r is 
the spatial distance between two points on the LSS sustained by an angle 9. Here C(r) is related with 
0f(?7o, k) via the angular power spectrum d in the following manner. 
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If we integrate by parts the right-hand side of Eq. (7) and use the Fourier-sine formula, we obtain a 
first-order differential equation for P(k), 



POO 

-f(k)k 2 P'(k) + [-f(k)f'(k)k + g 2 (k)] kP(k) =4tt / C* thin (r) sin krdr , 

Jo 



(9) 



where C thln (r) should be calculated using 9^ hin (J7o, k) in (8). Because of the oscillatory nature of the 
transfer functions, there appear zeros of f{k) where Eq. (9) becomes singular. However, the values of 
P{k) at these singularities are readily known as 



P(ki) 
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C thln (r) sin hrdr , 



(10) 



where fcj is the i-th zero point of f(k). Then the spectrum can be obtained easily and accurately despite 
the presence of the singularities, since we may solve it as a boundary value problem between singular 
points rather than an initial value problem. 

The thin LSS assumption, however, is not quite realistic and 0£ hm (?7o,fc) deviates from the exact one 
significantly. The thickness of the LSS must be taken into account for actual applications. In this paper, 
we present a new method, which is basically the same as the previous one but takes account of the 
thickness of the LSS, and hence can be applicable to realistic cases. 

We first improve the transfer functions in the formula (9) by including the effect of the thickness of 
the LSS. For this purpose, we use the following approximate expression instead of (5). 



0*™( m ,k) = (2£+l) 
+ (2£+l) 
+ (2£+l) 
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The first two terms of the above formula partially take account of the thickness of the LSS, which would 
be exact if the spherical Besscl functions would not oscillate but remain constant over the thickness of the 
LSS, which is the case for low wavenumber modes. Note that the ISW term, which was ignored in Paper I, 
is now included in the third term. It takes account of the early ISW effect, since the matter-radiation 
equality time is fairly close to the decoupling epoch. Although the width of the time interval during 
which the integrand of the third term is non-negligible is somewhat larger than the width of the LSS, 
the early ISW effect happens to be described by this approximation fairly accurately. In fact, Eq. (11) 
turns out to be a very good approximation for C(r) or Ce as is seen in Fig. 1, in which an angular 
power spectrum based on our approximate formula is compared with the fully time-integrated one. Note 
also that Eq. (II) is applicable to flat ACDM models as well, if we neglect the late ISW effect which is 
important only for low-multipolcs. In this case, we should not integrate the ISW term until 770, otherwise 
it causes serious error since the late ISW effect is incorrectly included. 
Defining new transfer functions f(k) and g(k) to express (11) as 



6^(770, k) = (21 + l)/(fc)$(0, k)j e (kd) + (21 + l)g(fc)$(0, k)j' e (kd), 



(12) 
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and repeating the same procedure as in Paper I, we arrive at the following differential equation. 



/ (k)k 2 P'(k) + -2f(k)f (k)k+f(k) kP(k) = A7r C app (r) sin krdr , 



(13) 



where the right-hand-side should now be calculated using (11) in (8). 

The above differential equation can be solved in the same way as (9) to obtain P(k). Unfortunately, 
however, if we used the observed correlation function obtained by, say, MAP satellite [3] in the right-hand- 
side of (13), we would reach an incorrect primordial spectrum P(k), because, although the approximate 
formula (11) reproduces the locations of peaks and troughs of the true spectrum quite well, the amplitude 
still deviates from the one obtained by full numerical calculation [13] which is to be compared with the 
observed one. 
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FIG. 1. The angular power spectra of CMB anisotropy in a spatially flat universe with h = 0.7, flj, = 0.03, 
f^CDM = 0.97 and scale-invariant adiabatic perturbation. Solid line is based on fully time-integration, while 
dashed line is the result of our approximate formula (11). The errors on large angular scales with small I are 
almost due to the approximation of the ISW term. On small scales (large I), the effect of thickness of the LSS 
becomes so significant that the approximation leads to large errors. However, the qualitative features such as 
peak locations are still reproduced relatively accurately. 



Fortunately, there exists a remedy for errors caused by the formula (11). To show this, let us consider 
the ratio of the exact to the approximated angular power spectrum, which we define as 

^~<cxact 

b e = -£mr- ( 14 ) 

Of course, this depends on the shape of the initial spectrum. However, the dependence turns out to 
be relatively weak as is illustrated in Fig. 2. That is, for a fixed set of the cosmological parameters, 
be is found to be rather insensitive to the variation of the initial spectrum. This suggests the following 
procedure of reconstruction. Consider we are given the observational data (= C^ hs ). Assuming the 
cosmological parameters are known, we calculate the correction factor be for a fiducial initial spectrum 
P(°\h). Namely, we calculate the exact angular spectrum C^p and the approximate one (7 app ( ) using 
Eq. (11) from P(°\k) and take their ratio which we denote by bf\ Then dividing Cf° s by bf*\ we 

estimate C app . Let us denote this estimated approximate angular spectrum by Cf . Then, we can 

insert it into Eq. (13) to reconstruct the initial spectrum. If be were rigorously invariant, this procedure 

would recover the true initial spectrum. Due to small errors caused by the non-invariance of be, however, 

the reconstructed spectrum, denoted by P^ (k), will not be exactly equal to the true primordial spectrum 

P(k). Nonetheless, we can expect P^(k) to be fairly close to the true P(k), or at least better than the 

fiducial spectrum P(°\k) which is a blind guess. Then, we can iterate this procedure to improve the 
accuracy. 

Schematically, this iterative procedure is described as 



pW(fc) - bP = -^-n - CV app(n+1) = ^ — > P(" +I )(fc), (15) 
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where the last step is the inversion procedure, and Cg is the approximate angular spectrum for 
using the formula (11). We repeat this procedure until we obtain a given degree of convergence. 
The validity of the above prescription and the rate of convergence depend on the degree of invariance 
of be. In Fig. 2, we plot Ce and be for a scale- invariant and a peaked spectrum. Although the difference 
of Ce between the two spectra is significant, the difference of be is small, namely, only a few percent. The 
invariance of be, of course, will not hold between entirely different initial spectra. Recent observations, 
however, suggest that the initial spectrum is almost scale-invariant on very large scales. We therefore 
assume that the scale-invariant spectrum is a good guess, apart from possible features like peaks and dips 
and/or some smooth variations of the power-law index over a range of one or two orders of k. 




III. RESULTS 



In order to test the validity of our method, we apply it to several shapes of initial spectra in flat CDM 
models. It may be noted that, although we focus on flat CDM models, we method is equally applicable 
to flat ACDM models. In fact, our models may be approximately regarded as representing flat ACDM 
models for the same value of Ocdm/i 2 [14]. This is because the angular spectrum depends on ficDM only 
through the combination Ocdm^- 2 , except for Ce at small values of £ to which the late ISW effect may 
contribute, but we do not use Ce at £ < 30 for the reconstruction. 

For initial spectra having distinct features, a spectrum obtained at the first round of inversion sometimes 
contains regions of negative values or has extraordinary sharp peaks or dips (with height or depth of more 
than a factor of ten relative to the continuum) due to errors in the estimation of be- Since a real spectrum 
must be positive definite, we make it a rule to cut off the negative part from the spectrum and smooth out 
sharp peaks and dips to mild ones, and use it for the next round of the inversion procedure. Although 
there is no justification for this prescription, wc find that these peculiar features disappear from the 
spectrum at the next round when they are spurious, while they show up again when they are real. We 
first assume the correct cosmological parameters are known. Then the iteration converges by a few times, 
although the rate of convergence depends on a spectrum. We find the converged spectrum agrees well with 
the original spectrum. In actual applications, we need to calculate the corresponding angular spectrum 
again to compare it with an observed spectrum C^ bs . 

We plot the result for a spectrum with a peak and a dip in Fig. 3(a). In this example, there appears a 
sharp dip at kd ~ 900 at the first round of inversion. Assuming it is spurious, we cut it off by hand and 
interpolate the spectrum smoothly to proceed the iteration. In the next round, it disappears completely. 
As seen from the right panel of the figure, the iteration converges remarkably fast. The small hump at 
kd ~ 400 disappears without trace, while the dip near kd ~ 800 grows deeper and approaches to the 
real shape. The bottom panel shows the relative error at the fourth round of iteration. Although the 
precision around the singularities of Eq. (13) is relatively bad, the error is still within 4%. We also plot 
the result for a spectrum which has smooth variations of the power-law index in Fig. 3(b). In this case, 
there appears no spurious peak or dip through the iteration so that we can easily recover the shape of 
the real spectrum without modification. 
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FIG. 3. Test of our inversion method for a primordial spectrum (a) with a peak and a dip superposed on a 
scale-invariant spectrum and (b) with smoothly changing power-law index, (a) The top left panel shows p' 1 ) (fc) 
which is obtained from the first round of iteration. The solid curve is the original P^'(k) which happens to have 
a very spiky dip at kd ~ 900. We assume such a feature is spurious and remove it by simply cutting that part off 
and interpolating from the both sides, which is plotted by the dotted curve. The top right panel shows P*- 2 ' (fc) 
and P( 4 \k) together with the original (true) spectrum P(k). The bottom panel shows P^(k)/P(k). The relative 
error is negligible except for regions near the singularities (marked by triangles) of the differential equation (13). 
Even near the singularities, the relative error is smaller than ~ 4%. (b) In this case, we do not need to modify 
P ( - 1 \k). The result reproduces the variation of the real spectrum well. 
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In the above, we assumed that the resolution of observations is limited to £ < £ c — 1500, which roughly 
corresponds to kd < k c d = 1500 in the primordial power spectrum. For our iterative procedure, however, 
we must assume the shape of the spectrum at scales k > k c , because the contribution of modes at k > k c 
to the angular spectrum at £ < £ c is non-negligible. It becomes small enough only for I < ^ c /2. Since 
there is no principle to determine the spectrum at k > k c , we simply extrapolate the spectrum at k < k c 
to k > k c with a power-law form. This induces some errors in the angular spectrum at £ > £ c /2- In the 
case of the example shown in Fig. 3, the original spectrum is assumed to be scale-invariant at k > k c with 
the amplitude extrapolated from that at k < k c . Hence we were able to recover the angular spectrum 
over the whole range up to £ = £ c . In an actual application of our method, however, we must decide not 
to use Cg hs at £ > £ c /2 as a criterion for the convergence. 



IV. CONSTRAINING COSMOLOGICAL PARAMETERS 

In the previous section, we showed that the initial spectrum can be determined with a high precision. 
However we had to assume that the cosmological parameters, which are usually estimated from the CMB 
data, are exactly known by some other observations from the beginning. Here we examine consequences 
of the use of incorrect cosmological parameters in our procedure. 

Naively, one might expect that our method is incapable of determining the cosmological parameters, 
since a different choice of cosmological parameters would simply give a spectrum that differs from the true 
spectrum. This is true in the rigorous sense. However, the very fact that the observed angular spectrum 
is used to determine the cosmological parameters indicates that there is not much freedom in varying 
cosmological parameters. In fact, if we use parameters that differ significantly from real values, very spiky 
features appear at the singularities of Eq. (13) in the reconstructed spectrum. Furthermore, depending 
on the direction of deviation from the real values, the reconstructed spectrum becomes negative in some 
regions of k, usually near the singularities, which never disappear by iteration. 

As an example, Fig. 4 shows the reconstructed spectra for several values of h and when Cf° s is 
given by a scale-invariant spectrum with h = 0.7 and Of, = 0.03 in the flat CDM universe. As noted 
before, we may regard this case to describe a flat ACDM models with f^cDM^ 2 = 0.7 as well. As for the 
h dependence, spiky peaks appear at the singularities for larger values of h, while spiky dips appear for 
smaller values of h and the dips become too deep to render the spectrum negative for h < 0.65. Therefore, 
the positivity condition of the spectrum severely constrains models with smaller values of h and larger 
values of fib, while models with larger values of h and smaller values of fib can be also constrained unless 
there is a good reason to believe that locations of the spikes and the singularities should coincide. 



V. CONCLUSION 

We have presented a method to solve the inversion problem of reconstructing the primordial curvature 
perturbation spectrum P{k) solely from the CMB angular power spectrum Ce- In Paper I, we developed 
an inversion method by deriving a first order differential equation for P(k) with its source term determined 
by Ce, but under the assumption of an infinitely thin LSS. In this paper, we improved it significantly by 
fully incorporating both the thickness of the LSS and the early ISW effect. This was made possible because 
of an empirical fact that the spectral ratio be of the exact, full-numerically calculated angular spectrum 
C| xact to an approximate angular spectrum C^ pp using an analytic formula for the CMB multipolcs (11) 
is relatively insensitive to variations of the shape of P(k). This fact allowed us to apply the inversion 
method developed in Paper I iteratively to reconstruct P{k). We have found the new method can recover 
the original spectrum accurately, to within relative errors of 4%, even in cases when there are distinct 
features like peaks and dips in the spectrum. 

Our method, however, has a possible drawback that the inversion procedure can be performed only 
when the values of cosmological parameters are known. Therefore, we have also studied the effect of using 
the cosmological parameters that are different from their true values. We found that a small deviation 
from the true values results in the appearance of peaks and dips at the locations of singularities of the 
differential equation for P(k), which can be easily judged as spurious. In particular, depending on the 
direction and the size of deviations, the spurious dips can become so deep that the positivity condition 
of P{k) is violated there. Thus, contrary to our concern about the inversion method mentioned above, 
it may be regarded as an advantage in the sense that the cosmological parameters can be constrained 
severely with no regard to the shape of the initial spectrum. This is in accordance with a widely accepted 
assertion that the cosmological parameters can be determined from the heights and locations of peaks in 
the observed CMB spectrum. What is new is that we have not only justified this assertion qualitatively 
but also provided a means to quantify the level of its validity or limitations without any assumptions on 
the form of the primordial spectrum. 
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Finally, let us point out a couple of issues to be explored in future studies. First, we should investigate 
how observational errors on Cg would be reflected to the reconstructed spectrum of P{k). Second, our 
method in the present form can apply only to spatially flat universes. To make it more general, an 
extension to cases of spatially curved universes is necessary. In this respect, we may note the following. 
Apart from the geometrical effect of curved space that changes an angle sustaining a fixed distance on 
the LSS, we may adopt the small-angle approximation which allows us to use various flat space formulas. 
Hence, an extension to non-flat universes seems feasible enough, if not straightforward. Another issue 
is about the CMB polarization. Our method assume that the CMB spectrum is dominated by the 
scalar-type curvature perturbations. However, it has been argued that the CMB spectrum may contain 
a non-negligible contribution from tensor perturbations [15]. To identify the tensor contribution to the 
CMB spectrum, the CMB polarization spectrum plays a crucial role [16]. Hence, to complete the inversion 
problem, it is necessary to develop a formalism that utilizes both the temperature and polarization spectra 
to reconstruct both the scalar and tensor perturbation spectra simultaneously. These issues are currently 
under study. 
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FIG. 4. The cosmological parameter dependence of the reconstructed spectrum. The real spectrum is 
scale-invariant with h — 0.7 and Qb = 0.03 in the flat CDM universe. We plot the reconstructed spectra for 
h = 0.65, 0.68, 0.7, 0.72 and 0.75, and Q b =0.02, 0.03 and 0.04. 
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